Atom-dimer p-wave resonance for fermionic mixtures with different masses 
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We show that, near a Feshbach resonance, a strong p-wave resonance is present at low energy 
in atom-dimer scattering for 6 Li- 40 K fermionic mixtures. This resonance is due to a virtual bound 
state, in the atom-dimer system, which is present at this low energy. When the mass ratio between 
the two fermionic elements is increased, this virtual bound state goes to a known real bound state 
which appears when the mass ratio reaches 8.17. This resonance should affect a number of physical 
properties. These include the equation of state of unbalanced mixtures at very low temperature but 
also the equation of state of balanced mixtures at moderate or high temperature. The frequency and 
the damping of collective modes should also provide a convenient way to evidence this resonance. 
Finally it should be possible to modify the effective mass of one the fermionic species by making 
use of an optical lattice. This would allow to study the strong dependence of the resonance as a 
function of the mass ratio of the two fermionic elements. 

PACS numbers: 03.75.Ss, 03.65.Nk, 34.50.Cx, 67.85. Lm 



One of the most remarkable features of the physics of 
ultracold gases is the simplicity of the effective interac- 
tion which is completely characterized by the scattering 
length between the various involved atoms. For example 
in fermionic gases involving only atoms in the two low- 
est hyperfine states of a single element, such as 6 Li or 
40 K, the interaction is fully described by the scattering 
length a between atoms belonging to these two different 
hyperfine states pQ. Remarkably this simplicity extends 
even to apparently more complex situations. Indeed on 
the BEC side of a Feshbach resonance, where a > 0, 
fermionic atoms in two different hyperfine states may 
form molecules, or dimers, with binding energy Eh- One 
might expect that the formation of these bosons could 
lead to more complex physics, expected for a molecular 
gas. Nevertheless the T = equation of state on this 
BEC side is very well described [I] in a wide domain 
of density by the Lee-Huang- Yang [5] equation of state 
for these composite bosons [3] , making merely use of the 
dimer-dimer scattering length 0,4 = 0.6 a between these 
bosons [H [5] . This is as if these bosons were elementary. 

In this paper we show that this simplicity is partially 
linked to the fact that the two atomic species have equal 
masses. It should be lost in the case of mixtures of atomic 
fermions corresponding to different elements, such as 
mixtures of 6 Li and 40 K toward which the interest has re- 
cently turned. Indeed new interesting physics is expected 
for large mass ratio [BTH2] . Quite recent theoretical work 
[l3[ IT4] has also focused on the case where the particle 
number is small and multiparticle bound states may oc- 
cur. Experimentally much progress has been made in 
exploring these mixtures [TMTT] and degeneracy has al- 
ready been reached. In particular a Feshbach resonance, 
which is fairly broad and reasonably stable, has been 
identified near 155G [18]. In the following we will as- 
sume that we are near such a Feshbach resonance and 



ignore any stability problem. We will follow the stan- 
dard use of refering to the two different atomic species as 
t and ! atoms, even if there is no relation between this 
notation and the physical spin of these atoms. We note 
a the scattering length between these atoms and assume 
that a zero-range interaction is a good enough approxi- 
mation so that the universal regime applies. Finally we 
ignore any trap effect. 

A first sign of more complex physics is the increase [19] 
(for a given a) of the dimer-dimer scattering length with 
the mass ratio m "f /m \,. However the most important 
point by far is the appearance of new bound states when 
this mass ratio increases. These bound states appear 
already at the level of the three-body problem and for 
the fermionic situation they have been studied in great 
details by Kartavtsev and Malykh [20], for two identi- 
cal fermions with mass m f and a different fermion with 
mass m J,. They are not present in the s-wave channel 
and appear only for angular momenta £ > 1. Among 
them are the well-known Efimov [21] states, which are 
quite remarkable for their spectrum structure with bind- 
ing energy going up to infinity when the range of the 
interaction goes to zero. For £ = 1 the Efimov states 
appear only fSOj [22] for a mass ratio m| fm \= 13.6. 
However there are already bound states for smaller mass 
ratios, which are not Efimov states. A first bound state 
[2Ti] appears for mass ratio m j" /m \.= 8.17. 

Actually we show in this paper that very important 
modifications in the fermion-dimer scattering properties 
arise for quite lower mass ratios. The basic reason is that, 
even before the mass ratio threshold for bound states is 
reached, virtual bound states are already present. For 
nonzero angular momenta, in particular £ = 1, the cen- 
trifugal barrier acts to inhibit their decay and provide 
them with fairly long lifetime. Hence they are not much 
different from real bound states. Moreover, since their 
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FIG. 1: Modulus of the £ = scattering amplitude j/o(fco)| 
as a function of fco, for mass ratios m^/ra^ — 1, 2, 4, 6.64 and 
8. The ratio mf/m; = 6.64 corresponds to the 6 Li - 40 K 
mixture. All these quantities are in reduced units, equivalent 
to set a — 1. 



energy is positive, they give rise to resonances, which are 
strong when they are located at low energy. We find 
that these effects are already important for the mass ra- 
tio m j" jm 1= 6.64 corresponding to the mixtures of 6 Li 
and 40 K, which is probably the most investigated exper- 
imentally. We note that these virtual bound states could 
play a very important role in the three-body decay and 
in the stability of these mixtures, since they correspond 
to physical situations where three atoms stay close to- 
gether for a fairly long time, increasing in this way the 
probability of three-body decay processes. We observe 
that the effect of the mass ratio can be studied experi- 
mentally, since the use of optical lattices allows a manip- 
ulation of the effective value of the atomic mass. These 
lattices could also be used to study and possibly remove 
instability problems. 

In this paper we are interested in mass ratios below 
the value 8.17, corresponding to the appearance of the 
first p-wave bound state. We consider only the s-wave 
and p-wave components of the full scattering amplitude, 
since the higher angular momentum components are ex- 
pected to give very small contributions to the total scat- 
tering. We show that the p-wave contribution becomes 
very strong, due to quasi-resonance, when the mass ratio 
approaches 8.17. 

This analysis is done quite conveniently by making use 
of the integral equation formulation of this dimer-fermion 
scattering problem first worked out by Skorniakov and 
Ter-Martirosian[23|. We use the full Skorniakov and Ter- 
Martirosian equations, corresponding to different fermion 
masses, for nonzero kinetic energy in the center of mass 
referential. 

We introduce the total mass of the dimer M = m t 
+m J,, its reduced mass /i = m^ra^j '(m-f + mj.) and its 
binding energy E^ = l/(2fia 2 ). We work in the center 
of mass referential for which the total momentum P = 
0. The corresponding total energy is E = — 1/(2/lmz 2 ) + 



FIG. 2: Modulus of the 1=1 scattering amplitude |/i(feo)| 
as a function of ko, for mass ratios m-f/mj, = 1,2,4 and 6.64. 
The ratio m^/m± = 6.64 corresponds to the 6 Li - 40 K mixture. 



fco/2/i-r where /it is the atom-dimer reduced mass fix = 
m^M/(m^ + M) = m^(m^+m^)/(2m^+m^), and ko the 
outgoing momentum of the f fermion. The masses enter 
only through the combinations R — (J./(J>t, R' = 2/i/mi. 

The unknown function 03 (k, k ) in the Skorniakov- 
Ter-Martirosian integral equation is, within a multiplica- 
tive constant, the T-matrix for this scattering problem. 
Its s-wave component a^ s {k,ko) is given by a3 S (fc, ko) = 
(l/47r) J dQ^aslk, k ), while its p-wave component 
reads a 3p (k,k ) = (l/47r) J dfl k Pi(k, k )a 3 (k, k ) where 
Pi(k, ko) = k • ko/fcfco is the I = 1 Legendre polyno- 
mial. The corresponding s-wave and p-wave scattering 
amplitudes are given by a,3 S (ko,ko) and a3 P (fco, ko) re- 
spectively. The integral equations satisfied by a 3s (fc,fc ) 
and a3 P (fc, k ) are obtained by projecting the general Sko- 
rniakov and Ter-Martirosian equation on the £ = and 
£ = 1 components. The resulting equations are quite 
similar. In order to make the equations simpler, we use 
systematically reduced variables: all the wavevectors are 
actually reduced (dimensionless) wavevectors, expressed 
in units of 1/a, and similarly a3 SjP (fc, ko) are reduced scat- 
tering amplitudes expressed in units of a. 

Introducing the notations: 



a{k, q) = 1 
(3{k,q) = R'kq 

we define the kernels: 



Rkl 



(1) 



a(k,q) + (3{k,q) 
a(k, q) — P{k, q) 



and 
K p (k,q) 



1 a(k,q) 
(3(k, q ) 2P{k,q) 



In 



a{k,q) +/3(k,q) 



a(k,q) - P{k,q) 



(2) 



(3) 



Then the equations for the s-wave and p-wave compo- 
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FIG. 3: Modulus of the ^ = 1 scattering amplitude |/i(fco)| as 
a function of fco, for mass ratios m^/m; = 6.64 (corresponding 
to the 6 Li - 40 K mixture), 7, 7.5, 7.75, 7.875 and 8. 



FIG. 4: Location of the maximum fcj of |/i(fco)| as a function 
of mass ratio m^/m^. 
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We consider first the results for the s-wave component. 
Since we are not interested in the phase shift, we restrict 
ourselves to the modulus |/o(fco)| = \^3s{ko, fco)| of the 
scattering amplitude. Here as well as in the following 
we investigate only negative values of the energy E in 
order to avoid consideration of inelastic processes, leading 
to the dimer breaking. This implies fcg < 1/R. The 
results for |/o(fco)| in terms of k are given in Fig[TJ for 
the mass ratios m^/m^ — 1, 2, 4, 8. We have also inserted 
the result for m^jm^ = 6.64 corresponding to the 6 Li - 
40 K mixture. 

The fco = value is equal to the scattering length, and 
increases from 1.179 corresponding to — to 2.099 
for = 8m^. This known fairly slow increase of the 
atom-dimer scattering length with the mass ratio is the 
dominant feature of Fig[TJ since otherwise all the results 
are of order unity for fco = 1 and keep decreasing for 
increasing fco. 

We turn now to the p-wave component results. We first 
give in Fig [2] the results for mass ratios m^/m^ = 1,2,4 
and 6.64, stopping at the value corresponding to the 6 Li - 
40 K mixture. For equal masses the p-wave component is 
always small, as it could be expected, and it can clearly 
be omitted when one deals with the atom-dimer scatter- 
ing properties. Hence the atom-dimer vertex can, to a 
large extent, be taken as a constant proportional to the 
atom-dimer scattering length. No complication is thus 
expected to arise from this side when many-body prop- 
erties will be investigated. The situation is qualitatively 
similar for m^/m^ ~ 2 and even 4: the p-wave compo- 
nent is fairly small and can be neglected compared to 



the s-wave contribution. However, we see that this p- 
wave component rises steeply when the mass ratio goes 
from 4 to 6.64. This last ratio is at the border of the 
resonant domain. 

This domain is now displayed in Fig. [3j where the re- 
sults for mass ratios m^/m^ = 6.64, 7., 7.5, 7.75, 7.875 
and 8. are plotted. We see that the hump, present for 
m t/ TO i = 6.64 around fco ~ 0.45, develops into a strong 
resonance at lower and lower energy when m^/m^ is in- 
creased. Although we have sampled mass ratios which 
are very close to each other, we see that the resonance 
grows very rapidly with increasing mass ratio. The phys- 
ical origin is clearly the development of a virtual bound 
state at positive energy. The lifetime of this state is 
directly related to the width of the resonance. Since 
the resonance peak gets quickly narrower, the lifetime 
grows very rapidly with increasing mass. The link with 
the bound state which appears [20 at zero energy for 
m t/ TO i = 8.17 is confirmed if we look at the position 
of the resonance peak as a function of the mass ratio. 
This is shown in Fig|4j We see that the position fcg of 
the resonance peak extrapolates to zero for a mass ratio 
which is in close vicinity of 8.17. More precisely when 
the resonance peak reaches fcg = 0, it will be infinitely 
sharp, corresponding to an infinite value for |/i(fco = 0)|. 
Hence in this case a^ v {k, 0) will be infinitely large, which 
implies that the homogeneous Eq.Q (i.e. without the 
term K p (k, 0) in the right-hand side) has a solution for 
zero energy fco = 0. This is just stating the well-known 
result that bound states are solutions of the homogeneous 
integral equation corresponding to Eq.Q. It is easy to 
find the lowest jra^ for which this homogeneous inte- 
gral equation has a solution. We find that this occurs for 
m^/mi = 8.172 in full agreement with Kartavtsev and 
Malykh gU]- 

Let us now come to the total scattering cross section 
cr(fco), which is likely to be the easiest direct physical 
quantity to measure experimentally. Since we neglect 
angular momenta higher than I = 1, we have for this 
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FIG. 5: Total scattering cross section (divided by 4-7r) Eq.|5| 
as a function of ho, for mass ratios rrif/mi — 1, 2, 4 and 6.64. 



FIG. 6: Total scattering cross section (divided by 47r) Eq.|5| 
as a function of fco, for mass ratios m-f/mj. = 1, 2, 4 and 6.64 
and 8. 



cross section: 

^=^(2f + l)|/,(fc )| 2 -|/o(fco)| 2 + 3|/i(A : o)| 2 (5) 

t 

The result is plotted in Figj5] for m^/m^ = 1,2,4 and 
6.64. We see that, while up to m^/m^ = 4, the cross 
section is almost featureless, it display a strong resonance 
for m^/m^ — 6.64 which corresponds to the experimental 
value for 6 Li - 40 K mixture. Naturally when the mass 
ratio is further increased, this resonance becomes even 
stronger. This is shown in Fig. [6] where we display also 
the case of m^/m^ = 8. The resulting resonance dwarfs 
the preceding one. 

This strong resonance in the atom-dimer scattering 
cross section for 6 Li - 40 K mixture should imply a number 
of experimentally observable consequences. The more 
natural way to evidence it would be to measure the colli- 
sion properties between a 40 K cloud and a 40 K - 6 Li dimer 
cloud, in a way analogous to the recent experiments on 
6 Li [24]. However this might not be so easy to perform. 
Measuring the equation of state of a mixture is on the 
other hand a fundamental question which should be eas- 
ier to answer. The resonance should affect in an im- 
portant way the low temperature equation of state of a 
mixture of 40 K cloud and a 40 K - 6 Li dimer, arising from 
an unbalanced mixture of 40 K and 6 Li atoms on the BEC 
side of the Feshbach resonance. This will not appear at 
the mean field level since the angular average will can- 
cel the effect because this is a p-wave resonance. But 
going to next order in 40 K, analogously to second order 
perturbation theory, there will be no cancellation. Hence 
one should see a strong dependence on the 40 K density. 
One expects also a strong density dependence to arise 
when the Fermi wavevector of the 40 K Fermi sea reaches 
the wavevector corresponding to the resonance, since for 
lower wavevectors the scattering is essentially negligible. 
One should see something similar to a threshold effect in 
density. We note that the resonance should also have an 
important effect even on a balanced mixture of 40 K and 



6 Li, on the BEC side, when the temperature is raised. 
Indeed the dimers will be partially broken by thermal 
excitation which provides a natural source of free 40 K 
atoms, and effects analogous to the ones arising in the 
unbalanced mixture. Similarly, going to high tempera- 
ture, we expect this resonance to have a marked effect 
on the virial coefficients since these are systematically 
related to the three-body problem [25] ■ Naturally the 
resonance is also expected to affect strongly the trans- 
port properties. Both the effects on the equation of state 
and on the transport properties should appear in the fre- 
quency and the damping of the collective modes [T] . 

Another way to understand the effect of this resonance 
on the Li - 40 K mixtures is to remark that the long lived 
virtual bound states responsible for it should behave in 
many respects in a way analogous to real bound states. 
Hence on the BEC side the physical description should 
involve not only free fermions and dimers, but also the 
existence of trimers. Obviously all the physical proper- 
ties should be affected by the presence of this additional 
fermion species. A related point is that the existence 
of these trimers should affect the dimer-dimer scattering 
properties. In other words we have only explored the sim- 
pler fermion-dimer scattering properties, but we expect 
that the dimer-dimer scattering amplitude will display 
related energy structure, and they should not be so com- 
plicated to explore with the methods of our preceding 
work [5j. 

Finally it is worth stressing that the use of optical 
lattices should provide a very convenient and powerful 
way to explore the resonance we have pointed out. In- 
deed the natural mass ratio of 6.64 between 40 K and 6 Li 
atoms happens to be just at the border of the strongly 
resonating mass ratio domain. A rather weak optical lat- 
tice could be used to slightly tune the effective mass of 
40 K or 6 Li (by making an appropriate choice of the light 
frequency, only one atomic species is affected), produc- 
ing strong modification of the resonance and hence of 
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the physical properties of the mixture. In this way one 
could go into the strongly resonating domain, or even 
reach the threshold for the appearance of real bound 
states. Or one could go in the other direction and ba- 
sically get rid of the resonance, which would allow to 
prove that it is responsible for specific physical prop- 
erties of the mixture. More specifically the 3D poten- 
tial F p t (r) = sE R [sin 2 (^ x) + sin 2 (K y) + sin 2 (iC z)] 
with E R = K^flm provides, within second order per- 
turbation theory, an effective m* given by m/m* — 
1 — ,s 2 /32. Increasing the 40 K effective mass to reach 
the bound state ratio 8, 17 requires s = 2.45 which cor- 
responds to a fairly weak optical potential. This simple 
picture of effective mass modification works only if the 
size of the involved objects (dimer, trimer) is large com- 
pared to the optical wavelength A. Taking A ~ 500 nm 
this should be a valid approximation in the vicinity where 
the scattering length becomes quite large. Even if correc- 
tions to this simple picture are necessary, the qualitative 
physical trends should remain valid. 

The authors wish to thank F. Alzetto for useful dis- 
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